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The effect of motion of heavy ions on the radiation characteristics of a point charge moving with 
uniform velocity along the direction of the external magnetic field in a plasma is investigated. For 
the case of stationary ions, two modes are found to be excited up to zero frequencies for values of the 
normalized strength of the external magnetic field above a certain minimum value. The result of 
inclusion of the motion of heavy ions is that the ordinary mode has a resonance at the so-called lower 
hybrid resonant frequency below which it is not excited and the extraordinary mode is excited up to 
zero frequencies even for very small values of the normalized strength of the external magnetic field. 
The power radiated in the ordinary mode in the neighborhood of the lower as well as the upper hybrid 
resonant frequencies is relatively large. The frequency and the angular spectrum of the emitted 
radiation as well as the direction of the Cerenkov rays are evaluated for some typical parameter values 
which include those usually obtained in the exosphere. 



1. Introduction 

A point charge moving with a uniform velocity in 
free-space does not radiate in view of the constraint 
that its velocity cannot exceed the phase velocity Co 
of the electromagnetic waves in the medium. On 
the other hand, a point charge moving in a material 
medium with a uniform velocity less than c can radi- 
ate provided its velocity is greater than the phase 
velocity of the electromagnetic waves in that medium. 
This radiation is the well-known Cerenkov radiation. 
An excellent treatment of the general field of Cerenkov 
radiation may be found in the book by Jelley [1958] 
and the review articles by Bolotovskii [1957, 1961]. 

A cold plasma in a magnetic field is equivalent to a 
dielectric, and Cerenkov radiation is possible in such 
a medium. However, in a cold plasma without ex- 
ternal magnetic field, the phase velocity of the elec- 
tromagnetic waves is always greater than Co with the 
result, Cerenkov radiation is not emitted in an iso- 
tropic plasma. But in an anisotropic plasma, Cerenkov 
radiation is possible. Kolomenskii [1956] and Sitenko 
and Kolomenskii [1956] have examined some aspects 
of the problem of radiation by a charged particle mov- 
ing with a uniform velocity along the direction of the 
external magnetic field in a plasma. In their papers, 
Kolomenskii and Sitenko have not systematically in- 
vestigated the frequency spectrum of the ordinary and 
the extraordinary modes nor have they given any 
consideration to the angular distribution of the radiated 
energy. Mckenzie [1963] has given a thorough treat- 
ment of the Cerenkov radiation in a magneto-ionic 
medium with emphasis on its application to the gen- 
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eration of low-frequency electromagnetic radiation in 
the exosphere by the passage of charged corpuscular 
streams. The present authors [Tuan and Seshadri, 
1963] have also given a treatment of the problem of 
radiation from a point charge moving with a uniform 
velocity along the direction of the external magnetic 
field in a plasma. It was found [Tuan and Seshadri, 
1963] that for certain parameters both the ordinary 
and the extraordinary modes were emitted up to zero 
frequencies. All the foregoing treatments of the 
Cerenkov radiation in a magneto-ionic medium are 
based on an idealization in which the motion of the 
heavy ions is neglected. It is reasonable to neglect 
the ion motion for frequencies sufficiently greater than 
the plasma and the gyromagnetic frequency of the ions. 
But for lower frequencies, it is not legitimate to ignore 
the ion motion. It is therefore of interest to study the 
effect of the motion of the heavy ions on the low- 
frequency spectrum of the radiation emitted by a 
uniformly moving charged particle in an anisotropic 
plasma. 

In this paper, the radiation characteristics of a point 
charge moving with uniform velocity along the direc- 
tion of the static magnetic field are investigated for 
the case of the plasma in which the motion of the heavy 
ions is taken into account. The method employed 
in this paper is different from those previously em- 
ployed [Mckenzie, 1963; Tuan and Seshadri, 1963]. 
In general, the emitted radiation is found to consist 
of two modes. The dispersion relations of these 
modes are analyzed in detail with a view to finding 
out the influence of the heavy ions on the low-frequency 
spectrum of the two modes. The frequency and the 
angular spectrum as well as the directions of the 
Cerenkov rays are examined. Numerical results are 
evaluated for some typical parameters of interest. 
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FIGURE 1. Geometry of the problem. 



In the frequency domain, the electric and the magnetic 
field vectors E(r, co) and H(r, go) are known to satisfy 
the following Maxwell's equations 



V X E(r, go) = i(OfXoH(F 9 to) (6) 

V X H(r, (o) =— iaj€ € • £(r, go) + zj z (7, go), (7) 

where p, and e are the vacuum permeability and di- 
electric constant. In (7), e is the familiar relative 
dyadic dielectric constant used in the magneto-ionic 
theory and its components are given by the following 
matrix: 



2. Statement of the Problem 

Consider a fully ionized, homogeneous plasma of 
infinite extent. For the sake of simplicity, the plasma 
is idealized to be a lossless and macroscopically neu- 
tral mixture of a gas of electrons and a single species 
of ions in which compressional effects due to finite 
temperature are neglected. An external magnetic 
field Bo is assumed to be uniformly impressed through- 
out the plasma in the z-direction, where p, </>, and z 
form a cylindrical coordinate system (fig. 1). It is 
desired to investigate the radiation characteristics 
of a point charge moving uniformly in the plasma me- 
dium along the direction of the external magnetic 
field. Let 
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represent a point charge g moving with a uniform 
velocity u along the z-axis from z = — » to z = °° such 
that it passes through the origin at t = 0. The current 
density arising from this uniformly moving charge is 
given by 



J (?, t)=zq u - — 8(z — ut\ 

llTp 



(2) 



where r represents the position vector of a point in 
the p, </>, and z-space. It is assumed that the source 
(2) is sufficiently weak so that the linearized plasma 
theory is applicable. 

Let E(T 9 t) and H(r, t) be, respectively, the electric 
and the magnetic field vectors. It is convenient to 
introduce the time Fourier transforms defined by 

/(r»=[_°° f(T,t)e^'dt (3) 

f(T,t)=±j X f(r,o>)e^»dw (4) 

to the source (2) and all the field quantities. The 
Fourier transform of the source (2) is obtained as 



J z (F,(o) = q ^e^l". 



(5) 
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where 
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at = l 
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(9a) 

(9b) 
(9c) 
(9d) 
(9e) 



In (9), a>pe is the plasma frequency, and a) ce is the gyro- 
magnetic frequency of the electrons. Also a> P i and 
cod are the corresponding quantities for the ions. It 
is found to be advantageous, though somewhat un- 
conventional, to employ the following normalized 
frequencies: 



fl=- 



Ope 



R= 2>2L. 

(Ope 



(10) 



With the above notation, (9a— c) becomes 



€1 = 1 



1 
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where 






(12) 



and m e and m\ are the masses of an electron and an 
ion respectively. 

In view of (5) and the geometry of the problem it is 
obvious that the field components are independent of 
(/> and depend on z only through the phase factor e iwzlu , 
which may be conveniently separated out as follows: 



f(r, cu)=/(p, o))e 



iioz/u 



(13) 



The arguments p, co of all the functions will be omitted 
for convenience. 

The projection of a vector in the p — </> plane will 
be denoted by the subscript t and z denotes the unit 
vector in the z-direction. The longitudinal and the 
transverse components of (6) and (7) are given by 



pdp 

i_d_ 

pdp 



{pEcfy) + (DPL{)H Z = 

(pHt) + iw€ {) €*E z = J z 
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(14a) 
(14b) 
(15a) 
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(15b) 



On cross multiplying (15a, b) by z, the expressions for 
zXH t and zXE f can be found and these when substi- 
tuted back in (15a, b) result in two simultaneous equa- 
tions in E ( and H t . The solution of these simultaneous 
equations yields the following expressions for Et and 
H t in terms of H z and E z : 



we 6 dp wue dp 

E^-iffMnft + M,*] 
dp 
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(16a) 
(16b) 
(16c) 

(16d) 



where 
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and Co — l/Vp, €o is the free-space velocity of electro- 
magnetic waves. 

The substitution of (16b and d) respectively in (14a 
and b) yields a pair of coupled wave equations, which 
are conveniently written in the matrix notation as 
follows: 



[M]V f[F] + [N][F] = [S] 



(18) 
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[S) = 
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(19c) 
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\M] = 



[N] = 



and N\ = wp, () , /V2 = weo6 ;i . (19d) 

The solution of (18) and the subsequent evaluation 
of all the physical quantities of interest are carried 
out by employing the methods first used by one of the 
authors [Seshadri, 1965b] in the treatment of the Ceren- 
kov radiation in a warm anisotropic plasma. 

3. Solution of the Coupled Wave Equations 

(18) 

The coupled wave equations (18) may be decoupled 
into two separate wave equations by a transformation 
to a new base. Let 



w]=mm- 



T n T22 



¥2 



(20) 



The substitution of (20) into (18) and the premulti- 
plication by [^-'[M] -1 yields 

V 2 m + [71-'[M]-'[/V][r]m = [r]- 1 [M]- 1 [S]. (21) 

The matrix [A/] _1 [7V] may be shown to possess two 
distinct eigenvalues k\ and k\, and hence, [T] may be 
chosen in such a way as to diagonalize [M] _1 [/V] with 
the following result: 



[V 2 + A:f,j^„=-S„7 



M 



2np 



for n = 1,2. (22) 
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With the help of (5), (13), and the right side of (21), it It was shown that 
is found that 



s, 


= -qo[T]^[M]- 1 





s 2 




1 



(23) 



The solution of (22) is easily seen to given by 
n (p) = ^S n H ( b m \k n p); /i=l, 2 and m = \ or 2. 



S(7, *) = £(r, *)X//(r, *). 



(25a) 



Therefore, the total power radiated in the radial direc- 
tion per unit area and per unit frequency interval is 
easily obtained with the help of (3), (4), (13), and (25a) 
to be given by 



S p (a>)=^Re E(p, a>)xR*(p, o>) 



= — Re [E+(p, co) - iH^ip, co)] 



From (16b, d), it follows that 
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(26) 



The substitution of (26) in (25b) and the use of (20) 
yields 



S P M = — Im 



dp dp _ 



[Ttf[M][T]* 



^2 



(27a) 



Since the matrix [A/] -1 [7V] is real, it follows that its 
eigenvalues k\ and k\ are either purely real or complex 



[T]-nM]- 1 [N][T] = 



k\ 
kh 



On premultiplying both sides of (27c) by 
[A] = [7f[M][T\ = 



an an 

CL21 CL22. 



(27c) 



(28) 



it is found that 



(24) [B] = 
The outward power radiated per unit area is given by 
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From (29) it follows that 

bii = ank\, b 22 = a 22 kl, 
bi2 = ai 2 kl, b 2 i = a2ik 2 v 



(29) 

(30a) 
(30b) 



It is known from (17) and (19) that [M] and [N] are 
symmetric, and therefore it follows that [A] and [B] 
are also symmetric. If the eigenvalues k\ and k\ are 
distinct, since 612 — b 2 \ and ai2 = #2i, it follows from 
(30b) that ai2 = a2i = bi 2 = b 2 i = 0. Consequently, 
it is found from (27b) and (28) that 



1 T d^i 

(25b) S p (co) = ^ Im I an** ^p + a 22 * 



*^2 

2 dp 



(31) 



Since the radiation is circularly symmetrical about the 
direction of motion of the charge, namely the z-axis, 
the total power radiated per unit path length and per 
unit frequency interval is obtained from (31) to be 
given by 



I(co) = 27TpS p (co) = p Im 



an** 



dp 



cW 2 



■****]■ ™ 



The total power radiated is seen from (32) to be the 
sum of the powers in the two modes separately. This 
arises due to the absence of cross terms in (32) and 
therefore the two wave functions denned by 1//1 and 



conjugates of each other. When the eigenvalues are ^2 are orthogonal. 

complex, the corresponding wave functions are seen 0n substituting (24) in (32) and making use of the 

from (24) to either increase or decrease exponentially, relation 



The exponentially grtrwing solution is ruled out from 
physical considerations and the exponentially decaying 
solution will not give rise to any radiation. The 
interest is therefore in the ranges of the parameter 
values for which the eigenvalues are real. For the 
case of real eigenvalues, [7] is real and hence (27a) 
becomes 



lmH^*{x)H^'(x) = (-\y 



7TX 



(33) 



it may be shown that 
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Sp(o))=7r-lm 
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For each mode, which corresponds to different values 
of n, m is chosen to be either 1 or 2 depending on 
whether the corresponding value of b nn is positive or 
negative. For a propagating mode, the value of A: 2 n is 
positive and with the above choice of m, I(co) given in 
(34) will be positive. Thus the fulfillment of the radi- 
ation condition, which requires a net outward flow of 
power at large distances from the source, is ensured. 

4. Dispersion Relations 

Since k\ and k\ are the two eigenvalues of [A/] -1 [TV], 
it is obvious that k\ 2 = V\la) 2 and k% 2 = V\lo) 2 are the 
eigenvalues of [</V] _1 [Af], where V is the velocity of 
phase propagation in the radial (p) direction. Using 
(19b, c), V 2 may be easily shown to satisfy the following 
quadratic equation: 



positive, real values of V 2 are seen to be obtained only 
in the following two cases: 
Case 1 B<0. 



this 



case, 



V(A/2) 2 -B |>|| 



V 4 -\-AV 2 ^-B = 9 



where 



and 
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B- 



2l Mn,M 2 2 



(MuM22-M 12 M2l). 
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(35) 



(36a) 



(36b) 



and hence V\ only is positive real, with the result that 
only the ordinary mode is excited in this range. 

I A\ 2 
Case 2 B > 0, A < and (-] > B. It can be easily 

shown that for this case V\ and V\ are botli positive 
real and hence both the ordinary and the extraordinary 
modes are excited. 

It is desired to find the ranges of the parameters 
ft and R (given by (10)) which correspond to the above 
two cases. In this way, the values of the parameters 
ft and R for which each of the two modes are generated 
can be determined. 

In most cases, the Alfven wave velocity V a — CQR\m 
is very much smaller than the free-space electro- 
magnetic wave velocity cq. Also the mass ratio m 
is very much smaller than unity. Consequently, in 
what follows, m and R 2 m are legitimately neglected 
in comparison with unity. 

The parameteric equations for A — oo, B = oo, A — 0, 
and B = may be easily shown to be given by the 
following curves in the ft 2 — J? 2 space: 



The dispersion equation (35) which specifies the phase 
velocity V for propagation in the direction perpen- 
dicular to that of the motion of the charge can also 
be derived in a much simpler way, as shown in the 
appendix. 

The use of (17a, b) and (19d) in (36a and b) yields, 
after some simplification, the following expressions 
for A and B: 



U: e 3 = 

U. i = ( 6l ±e 2 )£ 2 
L,: e,=0 

U 6i + €3 -j3 2 (€i€ 3 + €?-€$ = 



for A = °o and B = °° 
for A = oo and B = oo 
for£ = 
for .4=0. 



With the help of (11a, c), the curves L\ — La are shown 
to be given by the following parameteric equations 



and 


e 3 [l-2e,/3 2 + (e 2 -ei)/} 4 ] 


)1 

a (37a) 
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Q 2 =l+m 

[n 2 R* m+ { \ +m £ 2 ] 


= il 2 ft 2 (l-m) 2 


(39a) 
(39b) 
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"' (37b) 


Z, 3 :. 


Q 4 -n 2 (l+K 2 + m + /? 2 m 2 ) 




where 


C() 


(37c) 


Li^*. 


+ R 2 m(l + m + R 2 m) = 
2H 6 (1 — y3 2 ) - 2n 4 {(l + m) (1 - 2/3 2 ) 


(39c) 


The solutions of 

v\ = 

2 


(35) are given by 

A £ 


(38) 




+ R 2 (l + m 2 )(l -/3 2 )} + ft 2 [{2m 2 # 4 
+ (l + m)(l + m + m 2 )R 2 } 




-±±f^V(AI2)>-B 





where the subscripts 1 and 2 correspond to the posi- 
tive and the negative signs in (38) and denote ordinary 
and extraordinary, respectively. The factor €i/|ei| 
is introduced in front of the radical to ensure the con- 
tinuity of V2 across the line €i = 0. It is obvious that 
propagating waves are generated only in the ranges of 
parameters for which V 2 is positive real. From (38) 



-p 2 {R 2 m(l + m + R 2 m) + (I + m) 2 
+ (l + m)(l + m 2 )/? 2 
+ (1 + m + R 2 m) 2 }}-m 2 R\l + m) 
-r- P 2 R 2 m{\ +m) (1 + m + R 2 m) = 0. 



(39d) 



771 




Extraordinary, 
mode 2 




Figure 2a. Qualitative sketch of the regions of propagation of FIGURE 2b. Regions of propagation of the two modes . 
the two modes for a two component plasma. component plasma for /3 = u/c = 10~ 2 - 



In a similar manner, after some manipulations, it can 
be established that (AI2) 2 — B = is specified by the 
following parameteric equation: 

L 5 : R 2 [(l -P 2 ) (1 + m s )n 2 - m 2 (l + m) (l-(3 2 )R 2 

+ m{\ + m) 2 p 2 ] 2 = 4/3 2 (l - m 2 ) 2 (l + m- ft 2 )!! 4 . (39e) 

In figure 2a are sketched the curves L\—L$ which 
are seen to divide the Cl 2 — R 2 space into several re- 
gions inside which the signs of A, B, and (A/2) 2 — B 
do not change. The term A changes sign only on cross- 
ing the curves corresponding to A=0 9 °°; i.e., Li, L2, 
and L 4 . Similarly, B changes sign on crossing the 
curves corresponding to B = 0, °°; i.e., Li, Lz, and L 3 . 
Also, {A/2) 2 — B changes sign on crossing the curve L 5 . 
The signs of A, B, and (Aj2) 2 — B are determined for 
any one region, and then their signs for the other re- 
gions are obtained by inspection. It is then an easy 
matter to determine the regions of ll 2 — R 2 space cor- 
responding to the case 1 (B < 0, indicated by shading 
with vertical lines for the ordinary mode and with 
horizontal lines for the extraordinary mode) and the 
case 2 (B> 0, A < 0, and {A/2) 2 -B > 0, indicated by 
shading with horizontal and vertical lines). In figure 
2b, the curves L\— L5 are drawn to scale for the case 
yS = 10 -2 , and the regions of generation and propaga- 
tion of the ordinary and the extraordinary modes are 
indicated by shading with vertical and horizontal lines, 
respectively. 

In order to assess the effect of the motion of the heavy 
ions on the frequency ranges of excitation of the two 
modes, it is convenient to study briefly the limiting 
case of stationary ions, that is 771 = 0. For this case, 
the parametric equations of the curves L\ — L?> reduce 
to the following: 



L10: 2 = 1 



n 2 + 



1-/3 2 



= Vt 2 R 2 



(40a) 
(40b) 



L 3 p: O 2 = 0; H 2 = l+/? 2 
L 40 : ft 2 = 0; 



zi«= (o*-i) n 2 + 



£ 2 



1-? 



«-!) 



L 50 :R 2 = 



4/3 2 



(l-)3 2 ) : 



U-ft 2 ). 



(40c) 

(40d) 
(40e) 



In figure 3, the curves L\q — L^q are depicted. The 
regions of generation and propagation of the ordinary 
and the extraordinary modes are obtained as before 
and are indicated, as in the case of the two component 
plasma, by shading with vertical and horizontal lines, 
respectively. 

A comparison of figures 2a and 3 shows that the 
inclusion of motion of heavy ions considerably modifies 
the low-frequency end of the spectrum. The curve 
L3 is seen to have two branches, L3 and L'J, whose 
parametric equations are obtained to be given approxi- 
mately by n 2 = # 2 7rc/l + /? 2 -r-77i and H 2 = 1 + m + R 2 , 
respectively. It may be noted that Ll and U z represent 
the equations for the lower and the upper hybrid 
resonant frequencies, respectively [Stix, 1962]. It 
is seen that the motion of the heavy ions alters the 
radiation characteristics of a uniformly moving point 
charge in the neighborhood of the lower hybrid reso- 
nant frequency and lower. The lower hybrid resonant 
frequency is considerably larger than the ion gyro- 
magnetic frequency and is of the order of the geometric 
mean of the electron and the ion gyromagnetic fre- 
quencies. Since the lower hybrid resonant frequency 
can be in the ELF or the lower end of the VLF band, 
it becomes necessary to include the motion of the 
heavy ions in obtaining the low-frequency part of the 
Cerenkov radiation in a magneto-ionic medium. 

From figure 3, it is found that, in an electron plasma 
with stationary ions, the ordinary and the extraordinary 
modes are excited up to zero frequencies if the normal- 
ized strength of the external magnetic field is above a 
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certain minimum value which, however, is dependent 
on the particle speed. When the motion of the heavy 
ions is included, the ordinary mode is seen to have 
a resonance at the lower hybrid resonant frequency 
below which it is not emitted. The range of excitation 
and propagation of the extraordinary mode continues 
to extend up to zero frequencies no matter how small 
the normalized strength of the external magnetic field. 

Note that for the special value of R 2 = — 

m 1—p 2 

the extraordinary mode is not generated. As R 2 
moves away from the critical value, the frequency 
range of excitation of the extraordinary mode increases. 
Again for very small R 2 , the pass band of the extra- 
ordinary mode decreases as R 2 becomes smaller and 
finally vanishes for R 2 = 0. 

The normalized phase velocity Vi/co and the normal- 
ized wave number tlco/Vi for the ordinary mode and 
the corresponding quantities VzIcq and (lco/^2 for 
the extraordinary mode are evaluated numerically 
for four different values of R; namely, 1.0, 0.5, 
0.025, and 0.01, and are depicted in figures 4 to 
8. A perusal of the figures 4 to 8 reveals the fol- 
lowing interesting characteristics for the propaga- 
tion of the two modes in the radial direction. The 
ordinary mode has resonance, that is its phase velocity 
goes to zero, at the lower hybrid resonant frequency 
17 == V/? 2 m/l + R 2 + m and at the upper hybrid resonant 
frequency il — V 1 + m + R 2 . The phase velocity of 
the ordinary mode increases as the frequency is in- 
creased beyond the lower and decreased from the 
upper hybrid resonant frequencies, respectively. 
An examination of the frequency Cl versus wave num- 
ber ilco/V diagram shows that, in certain frequency 
ranges, the radial propagation is "backward" in char- 
acter in the sense that the phase and the group 
velocities are of opposite signs. In order for the total 
radiated power to be outwardly directed, it was found 
necessary to choose m=\ in certain frequency ranges 
and tfi = 2 in the other frequency ranges. It was 
found from the dispersion diagrams that the frequency 
ranges in which m=\ or 2 coincided exactly with the 
frequency ranges in which the phase and the group 
velocities are of the same or opposite signs, respec- 
tively. From (24) and the implied time dependence, 
it follows that m = \ corresponds to waves whose 
phase and group velocities are directed outwards from 
the source whereas m = 2 corresponds to a wave for 
which the group velocity is outwardly directed whereas 
the phase velocity is inwardly directed. 

It is found that the ordinary mode is a "forward" 
wave for H 2 > 1 -f m and is a "backward" wave for 
ll 2 < 1-h/n. The phase velocity of the ordinary mode 
decreases with the frequency for fl 2 >H-m but in- 
creases with the frequency for fl 2 < 1 + m. The extra- 
ordinary mode is found to be either "forward" or 
"backward" or partly "forward" and partly "back- 
ward" depending on the value of R. Its phase velocity 
is found to either increase {R = 0.025) or decrease 
(/\ = 0.01) with the frequency. It is interesting to 
note that in the magnetosphere, say at a distance of 10 



earth radii from the center of the earth, R is approxi- 
mately equal to 10" 2 and for that value of R, the phase 
velocity of the extraordinary mode decreases with the 
frequency. This feature can be easily seen to account 
for the "hydromagnetic whistlers" with a rising tone 
caused by the passage of the solar particles along the 
earth's magnetic field, in the magnetosphere. In a 
similar manner, the ordinary mode for il 2 < 1 + m with 
its phase velocity increasing with the frequency can 
account for the well-known atmospheric whistlers 
with decreasing tone. 
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FlGl kk 3. Regions of propagation of the two modes for a single 
component plasma. 




io"' n | 



FIGURE 4. Phase velocity and wave number of the ordinary and the 
extraordinary modes for R = 1 .0 and H. — 0.5. 
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and Fi, 2/co is given inJ38). The normalized frequency 
spectrum /i(o>) and h(aJ) are computed from (41a, b) 
for four different values of R; namely, /? = 1.0, 0.5, 
0.025, and 0.01. The normalized speed B=— is 

Co 

assumed to be equal to 10 -2 , which is the same as 
used in the calculations of phase velocity. It is of 
interest to note that the solar cloud speeds in the 
exosphere are of the order u= 10 _2 Co. The normalized 
frequency spectrum of the two modes are plotted in 
figures 9 to 17. A study of these figures shows that, 
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Figure 5. Phase velocity and wave number of the ordinary mode 
for R = 0.025 and SKI. 



5. Frequency and Angular Spectrum 

The frequency spectrum Ii(a>) of the ordinary mode 
and h(co) of the extraordinary mode are defined to be 
the total power radiated in the respective modes per 
unit path length and per unit frequency interval. It 
is of interest to examine the frequency spectrum of 
the two modes as a function of frequency for some 
typical values of R. With the help of (17b), (19d), 
(20), (23), (29), and (38), it is possible to show from 
(34), after some standard and straightforward manipu- 
lations, that the normalized frequency spectrum I\, 2(0)) 
is given by the following expressions: 
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Figure 6. Phase velocity and wave number of the extraordinary 
mode for R - 0.025 and R = 0.01 . 
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Figure 7. Phase velocity and wave number of the ordinary mode 
for R = 0.025 and Cl > 1. 
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Figure 8. Phase velocity and wave number of the ordinary mode for 
K = 0.01. 





FIGURE 10. Frequency spectrum of the ordinary mode for R = 1.0 
and 0< 1. 




Figure 9. 



Frequency spectrum of the extraordinary mode for 
R = 1.0. 
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FIGURE 11. Frequency spectrum of the ordinary mode for R- 
1.0 andil> I. 



in general, the power radiated is larger for the higher 
frequencies than for the smaller. It is seen that the 
frequency spectrum of the ordinary mode becomes 
infinite at the two hybrid resonant frequencies and on 
account of this, the total power radiated by a uniformly 
moving point charge becomes infinite. This difficulty 
can be overcome in the following manner. At the 
two hybrid resonant frequencies, the phase velocity 
in the radial direction goes to zero. The radial direc- 
tion is transverse to the direction of the external mag- 
netic field. In the neighborhood of the two hybrid 
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FIGURE 12. Frequency spectrum of the ordinary mode for R = 0.5 
andCK 1. 
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Figure 13. 



Frequency spectrum of the ordinary mode for R = 0.5 
and ft > 7 . 



resonant frequencies, since the phase velocity van- 
ishes, the magneto-ionic theory obviously fails. When 
the phase velocity becomes of the order of the thermal 
velocity of the charged particles, it is not legitimate to 
ignore the thermal motion of the particles. A treat- 
ment of the wave propagation transverse to the external 
magnetic field in a two component plasma has been 
given [Seshadri, 1965a] wherein the effect of the finite 
temperature of the electrons and the ions has been 
included. When finite temperature effects are in- 
cluded, it is found that the phase velocity instead of 
going to zero at the upper and the lower hybrid resonant 
frequencies, actually levels off to the value of the acous- 
tic velocity in the electron and the ion gas respec- 
tively. As a result of this, the frequency spectrum 
at upper hybrid resonant frequency, instead of becom- 
ing infinite, attains only a maximum value which 
depends on some inverse power of the acoustic veloc- 
ity u e in the electron gas. Similarly, the frequency 



spectrum at the lower hybrid resonant frequency 
attains only a maximum value which depends on some 
inverse power of the acoustic velocity u\ in the ion 
gas. Since u\ < u e , the frequency spectrum near the 
lower hybrid resonant frequency is likely to be higher. 
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Frequency spectrum of the ordinary mode for R = 0.025 
and 11 > 1. 




Figure 16. Frequency spectrum of the extraordinary mode for 
R = 0.01. 
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A systematic investigation of the Cerenkov radiation 
in a warm, anisotropic plasma has been carried out 
elsewhere by one of the authors [Seshadri, 1965b] for 
a single component plasma, and in that investigation 
a finite value has been obtained for the frequency 
spectrum near the upper hybrid resonant frequency, 
as anticipated above. 

The value of R over a Large part of the earth's exo- 
sphere is of the order IO -2 and for example, at 4 earth 
radii from the earth's center, /? = 2X10~ 2 , approxi- 
mately [Mckenzie, 1963]. It is seen from figure 2b 
that, for R between 2 X 10~ 2 to 6 X IO -2 , both the modes 
are emitted above the lower hybrid resonant fre- 
quency over a band which is largest for R = 2X IO -2 
and which continually diminishes as R is increased to 
the value R = 6X10~ 2 . It is seen that the power 
radiated in the ordinary mode is considerably larger 
than that in the extraordinary mode as can be seen 
from figure 14. In view of the presence of the lower 
hybrid resonant frequency, the power radiated in 
the neighborhood of that frequency is quite high, a 
fact which is not uncovered in the treatments wherein 
the ion motion is neglected. 

It is of interest to examine the polar diagram of 
the Cerenkov radiation. The normalized angular 
spectrum I\{0) and h(0) for the ordinary and the 
extraordinary modes are evidently given by the relation 



I„(d) =I„(il)l 2tt sin l} 



dO 



»=1, 2. 



(43) 



The angle between the direction of motion of the 
charge and the direction of the phase front is easily 
seen to be given by 



tan 6 H =p/(Vnlco) n = l, 2. 



(44) 



By making use of (38), (41), (43), and (44), the expres- 
sions for the normalized angular spectrum I\(0) and 
I'liO) of the modes can be obtained. The details of 
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FIGURE 17. Frequency spectrum of the ordinary mode for R = 0.01. Figure 18. Angular spectrum of the extraordinary mode for R = 1.0. 
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FIGURE 19. Angular spectrum of the ordinary mode for R = 
and Cl < 1. 



1.0 



these calculations as well as the final expressions, 
which are quite complicated, are omited here for the 
sake of brevity. The numerical results of the angular 
spectrum of the two modes are given in figures 18 
to 24 for the same four values of R as before and for 
/3=10~ 2 . The polar plot of the Cerenkov radiation 
shows interesting maxima and minima as a function 
of the angle made by the direction of the wave front 
with the direction of motion of charge. It is seen 
that the two hybrid resonant frequencies correspond 
to = 90°, as expected. Note from the figures 18 
to 24 that the angular spectrum has not been evaluated 
in certain cases up to very close to the hybrid resonant 
frequencies. 
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FIGURE 20. Angular spectrum of the ordinary mode for R = 1.0 
and Cl> 1. 
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6. Cerenkov Ray Direction 

In an anisotropic medium, the direction of the power 
flow or the Cerenkov ray is different from the direction 
of the wave normal. It can be shown [Stix, 1962] 
that the angle a between the directions of the Cerenkov 
ray and the wave normal is given by 



tan a n 



_ V n d /Co 
Co dO \V n 



n=l f 2 



(45) 



with the result the Cerenkov ray makes an angle 0-f- a 
with the direction of motion of charge. With the help 
of (38), (44), and (45), the Cerenkov ray direction speci- 
fied by + a corresponding to every wave normal direc- 
tion 6 can be computed. The numerical results are 
obtained for the same four different values of R as 
before and for the entire range of frequencies for which 
emission is possible. The numerical values as 



FIGURE 21. Angular spectrum of the ordinary mode for K — 0.5. 

depicted in figures 25 to 34 reveal two important fea- 
tures of the Cerenkov ray directions in an anisotropic 
plasma. In an isotropic medium, the Cerenkov ray 
directions make only acute angles with the direction 
of motion of the source whereas in an anisotropic 
medium such as the one considered in this paper, for 
appropriate parameter values, the Cerenkov ray makes 
obtuse angles with the direction of motion of the source, 
as may be seen from figures 26, 28 to 31, 33, and 34. 
The second important feature of the Cerenkov ray 
directions in an anisotropic plasma is that for ex- 
tremely low frequencies, the ray directions are con- 
fined to angles very close to the direction of the ex- 
ternal magnetic field, either in the same or in the 
opposite direction to the motion of the source. It is 
seen from figure 25 (extraordinary mode, 10~ 4 < £1 
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■FIGURE 22. Angular spectrum of the ordinary and the extraordinary 
modes for R = 0.025. 

<4.4X10~\ R = 1.0), figure 30 (extraordinary mode, 
10" 4 < ft < 10"\ R = 0.025), figure 31 (ordinary mode, 
6 X 10- 4 < a < 1 . 5 X 10-\ R = 0.025) and figure 33 (ex- 
traordinary mode, 10- 4 <H<2.3X10- 4 , «=0.01) 
that, for the frequency ranges and the appropriate 
mode type indicated within the brackets, the ray di- 
rections are confined to very small angles with the 
direction of the magnetic field. For larger values of 
R, the lower frequency Cerenkov rays are in the same 
direction as that of the motion of the source as seen 
from figures 25, 30, and 31. The atmospheric whis- 
tlers arising from the solar cloud streaming through 
the plasma of the exosphere in the direction of the 
earth's magnetic field are attributed to the ordinary 
mode (fig. 31) and have their ray directions in the 
direction of motion of the source. For smaller values 
of R, such as for example R = 0.0\, which are usually 
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FIGURE 23. Angular spectrum of the extraordinary mode for 
R = 0.01. 
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Figure 24. Angular spectrum of the ordinary mode for R = 0.01. 
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FIGURE 25. Ray direction of the extraordinary mode for R = 1.0. 



779 



180 

I5tf- 

a + 9 




R = 1.0 




J I I L 



10° 20° 30° 40° 50° 60° 70° 80° 90° 



FIGURE 26. Ray direction of the ordinary mode for R = 1.0 and 

n<i. 
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Figure 29. Ray direction of the ordinary mode for R = 0.5 and 

n> 1. 
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Figure 28. Ray direction of the ordinary mode for R = 0.5 and 
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Figure 30. Ray direction of the extraordinary mode for R = 0.025. 
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Fk;ure 32. Ray direction of the ordinary mode for R — 0.025 anil 
il> 1. 

obtained at distances greater than about five earth's 
radii from the center of the earth, the extremely low- 
frequency Cerenkov rays are confined to small angles 
about the direction of the magnetic field but lie in a 
direction opposite to that of the motion of the source. 
The sub ELF or the hydromagnetic emissions [Tepley, 
1961] are attributed to the extraordinary mode (fig. 33) 
and have their directions behind the motion of the 
source. It is seen from figures 25 to 34 that, in certain 
cases, the ray directions change very rapidly with the 
wave normal direction and the physical explanation 
for this behavior is not clear. 

To facilitate easy reference, an index to all the 
figures is provided in table 1 wherein for each value 
of R, the frequency ranges of excitation, the mode 
number, the wave type and the numbers of figures 
showing the dispersion, the frequency, and the angular 
spectrum and the ray directions are included. 
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Figure 33. Ray direction of the extraordinary mode for R = 0.01 . 
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FIGURE 34. Ray direction of the ordinary mode for K = 0.01. 
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7. Concluding Remarks 

The characteristics of some of the electromagnetic 
emissions which are detected both on the earth and 
by satellite borne receivers suggest Cerenkov type 
radiation generated by charged particles moving at 
high speed along the magnetic field lines in the 
magnetosphere as a possible mechanism. With the 
purpose of understanding these emissions, consid- 
erable investigations of the Cerenkov radiation in a 
uniform magneto-ionic medium have been carried out. 
In all these treatments, the ion motion has been 
omitted. Whereas the neglect of ion motion may be 
legitimate for the higher end of the VLF band, it is 
only proper to take into account the motion of heavy 
ions for the evaluation of the electromagnetic emissions 
in the lower end of the VLF and ELF bands. Re- 
cently, there has been interest in the interpretation 
of the observed emissions at hydromagnetic frequen- 
cies in which range the ion motion certainly plays a 
dominant role and cannot therefore be neglected. 

In this paper, the restriction on the motion of the 
ions is removed and the theory of Cerenkov radiation 
in a magneto-ionic medium is extended to apply even 
to ELF and hydromagnetic frequencies by the inclu- 
sion of the motion of the heavy ions. Not only has the 
analysis of the Cerenkov radiation in a two component 
magnetoplasma been carried out using a different and 
somewhat simpler method but some of the special 
features, such as, for example, the existence of the 
"backward" wave character of the radial propagation, 
which have been overlooked previously in the liter- 
ature, are emphasized in this treatment. 

In an anisotropic medium, the direction of power flow 
and the wave normal direction do not coincide. The 
radiation condition requires that the net power flow 
be directed outward from the source. In certain 
frequency ranges, it is possible for the wave normal 
direction to be inwardly pointed, giving rise to back- 
ward waves. It was found necessary while evaluating 
the frequency spectrum to require the waves to have 
inwardly traveling phase fronts for certain frequency 
ranges to obtain outwardly traveling power. The 
existence of the backward waves was further confirmed 
by the examination of the dispersion diagrams, which 
yielded oppositely directed phase and group velocities 
in the frequency ranges appropriate to the backward 
wave regions. 

Extensive numerical results on the frequency spec- 
trum, the angular spectrum, and the Cerenkov ray 
directions are presented for parameter values appro- 
priate to the magnetosphere and for a particle speed 
of the order of the estimated solar cloud speeds. The 
power radiated by a single point charge in uniform 
motion is so small that it cannot account for the ob- 
served power levels of low frequency emissions. A 
realistic estimate of the number density of particles 
whose coherent emission alone can give the proper 
order of magnitude for the intensity of the low- 
frequency emissions will require a detailed examina- 
tion of the propagation mechanism and the inclusion 
of the effect of the ionosphere boundary, in addition 



to the mechanism of excitation which alone is consid- 
ered in this treatment. 
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8. Appendix 

The dispersion relations for a plane wave propa- 
gating at an angle 6 with respect to the direction of 
the external magnetic field are given by the following 
Appleton-Hartree formula [Stix, 1962] 
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^o = € 3 (ef-e|) (A.2) 

At =- [€ + eie 3 - (e - 6ie 8 ) cos 2 0] (A.3) 



A 2 = €i — (ei— e 3 )cos 2 



(A.4) 



and v P h is the phase velocity of the plane wave in the 
magneto-ionic medium. If the objective is merely 
to determine the ranges of the parameters Cl and R 
in which the propagating waves are excited by a point 
charge moving with uniform velocity u in the direction 
of the external magnetic field, then it is possible and 
also desirable to start with (A.l) and obtain the disper- 
sion relation (35) specifying the phase velocity v for 
propagation in the direction perpendicular to the direc- 
tion of motion of the charge. It is well known that a 
charge moving with uniform velocity u through a 
medium of refractive index cjv p h will emit Cerenkov 
radiation at an angle with respect to its direction of 
motion if the condition for coherent radiation, namely, 



u cos 6 = v P h 



(A.5) 



is satisfied. Also, it is easily seen that the phase 
velocity v normal to the direction of motion of charge 
is given by 

v sin 6 = v P h. (A.6) 

With the help of (A.5) and (A.6), it may be shown that 

v 2 ph = v 2 /(l + v 2 lu 2 ). (A.7) 

The elimination of cos 6 and v p n from (A.l) — (A.4) 
with the help of (A.5) and (A.7) yields (35), after some 
simplification. 
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